Explanation for the isotropy of the Dirac cone in graphene 
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It is shown that in the absence ol spontaneous symmetry breaking the Dirac cones in the system ol 
interacting electrons on honeycomb lattice are isotropic at low energies. The effect is due to the Z3 
subgroup of the Dz symmetry group of the dispersion relation of Dirac quasiparticles. Consequences 
of the violations of the or the sublattice symmetry are discussed. 



The crucial property of the single layer of graphite, i. e. 
graphene, is that it has well-defined, Dirac-type, massless 
quasiparticles in its low-energy spectrum. It is commonly 
assumed that the Dirac cone is isotropic, with the Fermi 
velocity independent of direction. This is, for example, 
what one straightforwardly finds in simple tight-binding 
model of non-interacting electrons which is used as 
an input in the theories that include the Coulomb inter- 
action. Since all the components of the electron-electron 
interaction, including the Coulomb's unscreened ~ 1/r 
tail, are technically irrelevant perturbations at the non- 
interacting fixed point [2], it seems natural then to as- 
sume that the Dirac cone stays isotropic in presence of 
interactions as well. 

And indeed, the above expectation seems to be fulfilled 
in experiment Q. The above simple reasoning for the 
isotropy of the Dirac cone, however, is faulty. The irrele- 
vance of interactions is only a long distance property, and 
it is not obvious why the lattice effects, which need to be 
inserted by hand into the low energy theory, and would 
certainly renormalize the Fermi velocity, should not also 
break the cone's isotropy. The honeycomb lattice, after 
all, has only Cqv symmetry. Clearly, if an anisotropic 
cone would be assumed in the low-energy theory, for ex- 
ample, the interactions, being irrelevant, would not alter 
that hypothetical feature either. 

The issue is further highlighted by a recent study [4] 
in which the effects of weak interactions in a Hubbard- 
type model on honeycomb lattice are determined in a 
more rigorous form as a convergent series. The authors 
conclude that whereas the Fermi velocity does become 
renormalized by the short-range interactions, it does so 
in an isotropic fashion. This result comes out in a rather 
non-trivial analysis of the terms in the perturbation se- 
ries. The utter simplicity of the conclusion, and its obvi- 
ous relevance for the physics of graphene, however, beg 
for an elementary explanation. Such an explanation is 
offered here by showing that the single-particle excita- 
tion spectrum at low energy has the exact D3 symmetry, 
which suffices for the Dirac cone to be isotropic asymp- 
totically. The anisotropy may be induced by the reduc- 
tion in Hamiltonian's symmetry. The argument is based 
only on the existence of the well-defined quasiparticle ex- 
citations around Dirac points, and on the absence of the 
spontaneous symmetry breaking by the ground state. 

Let as assume that there exists a Dirac point in the 
Brillouin zone at a wavevector K near which quasiparti- 
cles are well defined in the sense of Landau. Cg^ sym- 



metry of the interacting Hamiltonian on the honeycomb 
lattice implies then that the quasiparticle spectrum e{k) 
is such that 



and 



eik + R) = eik), 



e{U^k) = e(fc). 



(1) 



(2) 



where R is any of the two reciprocal lattice vectors of the 
triangular Bravais lattice, and Ucf, is a rotation by the 
integer multiple of 7r/3. (We will subscribe here to the 
notation in [5].) This in particular dictates that there 
must be six Dirac points, related to each other by the 
above discrete rotations. 

Consider a point ki in the vicinity of the chosen Dirac 
point at K. There are two other symmetry-related points 
near the same Dirac point: 



^2 = t^27r/3^1 + Rl, 



(3) 



(4) 



so that obviously e{ki) = e(fc2) = e(fc3). These three 
points are related to each other by the rotations by 27r/3 
around the Dirac point. The Dirac cone therefore inherits 
the Z3 subgroup of the lattice symmetry. Furthermore, 
since for any wavevector k we also have the reflection 
symmetries. 



e{kx, ky) — e( k^, ky) — e{kx, ky), 



(5) 



Ky, then the 



it follows that, in the reference frame in which the x- 
axis is chosen along one of the reciprocal vectors 5] , one 
of these applies also to the vicinity of the Dirac point. 
More precisely, if we define the momenta relative to the 
Dirac point as Qx = kx — Kx and qy = ky 
refiection symmetry together with Eq. (1) implies 

e{qx + Kx,qy + Ky) = e{-qx + Kx,qy + Ky). (6) 

The full symmetry of the Dirac cone is therefore 

D3 = Z2{re flection) x Z^{rotations), (7) 

i. e. the symmetry of the equilateral triangle. 

It is convenient to define the complex wavevector mea- 
sured from the Dirac point: 



(8) 
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q = qx-iqy, (9) 

and place the origin of the coordinate frame in the recip- 
rocal space at the Dirac point, so that e{q + K) — > e((j). 

We will assume the particle-hole symmetry of the spec- 
trum. In graphene, the particle-hole symmetry is weakly 
violated by the very small next-nearest-neighbor hop- 
ping. At sufficiently low-energy, however, strictly speak- 
ing the particle-hole symmetry always emerges, after the 
chemical potential has been tuned to the Dirac point. 
Since in this note we are interested only in the asymptotic 
form of the energy-momentum relation, possible lack of 
particle-hole symmetry at higher energies does not affect 
our discussion in an essential way. 

Particle-hole symmetry of the Dirac spectrum implies 
that one should consider the square of the energy as a 
power series 

e^((j, q) — a\q + h\q + a2q^ + b2qq + C2q^ + (10) 
asq^ + bzq^q + c^qq^ + d^q^ + h.o.t., 

where Xi, x — a,b,c,d are constant coefficients. The 
symmetry of the spectrum implies that 

e{e^-'^/\,e-^-'^^'q)=e(q,q), (11) 

for n — 1,2, so that in the above expansion 

ai = 6i = a2 = C2 = 63 = C3 = 0. (12) 

Together with the reflection part of D3, which demands 
that 

e{q,q) ^ e{-q,-q), (13) 

this implies that 

e^{q, q) = b2iqq) + 03(9^ - f) + Ci{qqf + 0{q^q, q(f ), 

(14) 



where we included the single allowed quartic term as well. 
The reality of the expression requires the coefficients 62 
and C4 to be real, and 03 to be purely imaginary. 

To the leading order the spectrum near the Dirac point 
is therefore isotropic, and one can identify the coefficient 
62 as the (squared) Fermi velocity. 

It is also interesting to consider the situation where 
the Hamiltonian is deformed so that the spectrum near 
the Dirac point remains only Z2 (reflection) symmetric. 
This, for example, would arise if the nearest-neighbor 
hopping in one direction is different from the other two 
in the tight-binding model The energy then has the 
form 

e^(g, g) = ai{q ~ q) + 02(9^ + g^) + b2qq + h.o.t. (15) 

so that, besides the shift of the location of the Dirac point 
(due to a finite coefficient ai), the effect is the anisotropy 
of the Dirac cone (due to a finite coefficient 02). 

Finally, breaking of the sublattice or the time-reversal 
symmetry Q would only result in the addition of a fi- 
nite term to the final form of the Taylor expansion of 
the energy in Eq. (14), since the group remains the 
symmetry of the energy dispersion in the presence of the 
mass-gap. 

To conclude, whenever the Dirac quasiparticles appear 
in the spectrum of the interacting Hamiltonian that re- 
spects the symmetry of the honeycomb lattice, they in- 
herit the Dt, subgroup of the lattice symmetry group. 
The Z3 subgroup of the D3, when not broken, then forces 
the Fermi velocity to be isotropic. The anisotropy ap- 
pears only in the sub-leading term in the dispersion re- 
lation. 
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